Low- Temperature Spin Dynamics of Doped Manganites: 
roles of 'M.n-t2g and Cg and 0-2p states 

Priya Mahadevan^, I. V. Solovyev^'^ and K. Terakura^''* 
^ J RC AT -Angstrom Technology Partnership, 1-1-4 Higashi, Tsukuba, Ibaraki 305-0046, Japan 
'^Institute of Metal Physics, Russian Academy of Sciences, Ekaterinburg GSP-170, Russia 

^JRCAT-NAIR, 1-1-4 Higashi, Tsukuba, Ibaraki 305-8562, Japan 
^Institute of Industrial Science, University of Tokyo, 1-22-1 Roppongi, Minato-ku, Tokyo 

106-8558, Japan 
(February 1, 2008) 

Abstract 

The low-temperature spin dynamics of doped manganites have been analyzed 
within a tight-binding model, the parameters of which are estimated by map- 
ping the results of ab initio density functional calculations onto the model. 
This approach is found to provide a good description of the spin dynamics 
of the doped manganites, observed earlier within the ab initio calculations. 
Our analysis not only provides some insight into the roles of the Cg and the 
t2g states but also indicates that the oxygen p states play an important role 
in the spin dynamics. This may cast doubt on the adaptability of the con- 
ventional model Hamiltonian approaches to the analysis of spin dynamics of 
doped manganites. 
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There has been resurgence of interest in transition metal oxides with the perovskite 
structure owing to their wide range of electronic and magnetic properties. Among them, the 
hole doped manganites have been occupying a special position: they exhibit dramatic 
phenomena like colossal magnetoresistance and are being intensively studied with prospect 
for technological applications. LaMnOa, the parent material of the manganites, is an anti- 
ferromagnetic insulator. Upon sufficient doping {x ~ 0.15) with divalent ions (such as Sr, 
Ca), the system is driven metallic. The holes are allowed to move only if adjacent spins 
are parallel, which results in a dramatic increase in the conductivity when the spins order 
ferromagnetically, an effect which can be induced by applying a magnetic field or by lowering 
the temperature below the Curie temperature, T^. Thus, the carrier mobility is intimately 
related to the underlying magnetic state of the system, and there have been considerable 
efforts in recent times to identify the interactions that control the magnetoresistive prop- 
erties. An approach in this direction has been to analyze the spin dynamics of the doped 
manganites. 

Early experiments on Lao.yPbo.sMnOa 0] indicated that the spin- wave dispersion w(q) in 
the doped manganites could be interpreted in terms of a conventional Heisenberg ferromag- 
net with only the nearest-neighbor exchange coupling. This behavior is consistent with the 
double-exchange limit of the one-band ferromagnetic Kondo lattice model P], implying that 
conduction eg electrons move in a tight-binding band with one orbital per site and interact 
with localized t2g spins via the large intraatomic exchange Jh- However, more recent exper- 
iments on other doped manganites have found strong deviation of uj{q) from the simple 
cosine-like behavior expected from a nearest-neighbor Heisenberg model. Farther neighbor 
interactions in addition to the nearest neighbor one had to be taken into account to repro- 
duce softening of the dispersions for the wave vector q approaching the zone boundary. The 
one-band models could not explain the observed zone boundary softening, even qualitatively 
It was then suggested that additional degrees of freedom, probably the lattice degrees 
of freedom, may play an important role in the spin dynamics of these materials. Recently 
P], it was shown that the softening at the zone boundary has a purely electronic origin, and 
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could be explained within the framework of ab initio density functional band calculations, 
in the local-spin-density approximation (LSDA). 

The previous work also carried out a perturbative analysis of the exchange interaction 
strengths within a tight-binding model considering the double degeneracy of eg levels on the 
Mn site. It was argued that the degeneracy of Cg orbitals plays important roles, and simply 
by taking into account the proper structure of the kinetic hopping between nearest-neighbor 
Cg levels one can, to a large extent, understand the behavior of two strongest interactions, 
Ji and J4, in the half-metallic regime. Here, Jk corresponds to the exchange interaction 
between the A;-th neighbor atoms 0, as defined later by Eq. ^. Furthermore, it was pointed 
out that a realistic model including the oxygen p-orbitals and the Mn t2g orbitals could indeed 
modify the quantitative aspects of the results. While the fully-filled majority spin t2g orbitals 
could contribute an antiferromagnetic superexchange component to Ji, the partially filled 
minority spin t2g orbitals could, as suggested by the ab initio band structure calculations, 
contribute a ferromagnetic double exchange component. Further J2 was found to increase 
quite strongly in the e^-only model, unlike the weak dependence seen in the ab initio results. 
The previous work suggested that the oxygen p bands could modify J2 considerably as the 
Mn 3d-0 2p energy separation is comparable with the exchange splitting of the majority and 
minority spin orbitals. In the light of these observations, we have attempted to make further 
quantitative analysis to understand the origin of the observed zone-boundary softening. This 
has been done by mapping the results of the ab initio band structure calculations onto a 
tight-binding model which gives us flexibility of constructing simpler models and analyzing 
the contributions to the observed softening. 

The band structure for hypothetical cubic ferromagnetic LaMnOa with the lattice pa- 
rameter of 3.934 A, calculated within the linear-mufhntin orbital method with the atomic 
sphere approximation (LMTO-ASA), was mapped onto a nearest neighbor tight binding 
model p] which had been found to give a good description of the electronic structure of the 
transition metal oxides of the form LaMOs, where M=Ti-Ni. The tight-binding Hamilto- 
nian consists of the bare energies of the transition metal d (e^) and the oxygen p (ep) states 
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and hopping interactions between the orbitals on neighboring atoms. The nearest neighbor 
hopping interactions were expressed in terms of the four Slater Koster parameters, namely 
ppcr, ppn, pda and pdrr. Note that no direct d-d hopping was taken into account. While p-d 
covalency effects lift the degeneracy of the d orbitals, an additional interaction sda between 
the transition metal d and oxygen 2s orbitals was required to lift the degeneracy at the F 
point The energy of the oxygen 2s level was fixed at —20 eV. In order to obtain the 
magnetic ground state within the single particle tight binding model, we have introduced 
an extra parameter, (epoi) which is the bare energy difference between the up and down spin 
d electrons at the same site. An additional splitting, {e'p^i) was introduced between the up 
and down spin d orbitals of Cg symmetry The parameters entering the tight-binding 
Hamiltonian were determined by the least square fitting of the energies obtained from tight- 
binding calculations at several fc-points to those obtained from the LMTO calculations. It 
should be noted that the deep lying oxygen 2s bands were not involved in the fitting. The 
extracted parameters are sda = —1.57 eV, ppa= 0.91 eV, pp-K = —0.23 eV, pda = —2.02 eV, 
pdn= 1.0 eV, — £^=0.48 eV, epoi=S.2 eV, ep^i=0.3 eV being consistent with the earlier 
estimate for the system 

The frozen spin spiral approximation [|n|], where the orientation of the magnetic moment 
at each atomic site is spirally modulated by the wave vector q, was used to calculate the 
exchange interaction Jq defined by 

Jq = ^ Jfc exp[iq ■ Rfc] (1) 

k 

where Rj is the position vector of z-th Mn atom. Jj. is the k-th neighbor exchange interaction 
appearing in the Heisenberg Hamiltonian given by 

E[{ei}] = -^kGi ■ ei+k (2) 

ik 

with Gj denoting the direction of the magnetic moment at the site i. By using the local 
force theorem [|T2|, the changes in the single particle energy could be related to the exchange 
interaction by mapping onto the Heisenberg model as defined above. A rigid band picture 
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was adopted to simulate the doping effects [|T3[. Simplified models were constructed to 
elucidate the mechanism of zone-boundary softening of the spin wave, and comparison was 
made with the results from the LMTO calculations whenever possible to ensure that the 
present result is not an artifact of a particular parameter set. 

In Fig. la we show the LMTO results for the spin dispersion a;(q; x) along the symmetry 
directions FX, XM and MR calculated for several doping values x. In the small q region, the 
spin excitations have a weak dependence on doping as is evident from the result along the FX 
direction. However, sufficiently away from the F point, the results become a strong function 
of the concentration x. Considering the result for x=OA along the FX direction, we see that 
the spin dispersion is almost fiat from midway to the zone boundary. The experimental 
result for Pro.esSro.srMnOs along FX shows very similar behavior. The results of the 
tight-binding model, calculated by using the parameters extracted by fitting the ab initio 
band structure are shown in Fig. lb. This model calculation is called model A in order to 
distinguish from other models discussed later. Model A is seen to provide a good description 
of the energetics of the spin dynamics observed within the ab initio approach. The results 
in Fig. 1 suggest that the difference in energy between the ferromagnetic ground state and 
the various commensurate antiferromagnetic (AF) spin configurations such as those defined 
at the X (A-type AF), M (C-type AF) and R (G-type AF) points decreases with doping. 

The Fourier transform of Jq gives us the real space exchange integrals Jj as shown 
in Fig. 2. Dominant interactions are all confined within the linear -Mn-O-Mn-. . . chains 
parallel to < 001 > as was pointed out already 0. They are Jj with i=l, 4 and 8. J2 is 
the interaction for the pairs along < 110 > and takes relatively small values partly because 
of the cancelation between the contribution from the Mn d bands and that from O p states. 
These results are consistent with the analysis of the experimental results 0] which required 
finite J4 and Jg to be included in the Heisenberg Hamiltonian in order to reproduce the 
experimentally observed spin wave dispersions. In order to understand the behavior of the 
spin- wave dispersion ^(q; x) in terms of Jj, the following expressions for q || x will be useful. 
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1 3 

hu{q^) ~ 2[(Ji + 4J2) sin^ -q^a + J4 sin^ q^a + Jgsin^ -^qxo], (3) 

where a is the lattice constant of the cubic unit cell. For q^a -C 1 the above expression 
reduces to 

hu{q,) - ^[^1 + + 4J4 + 9J8](gxa)' (4) 

which helps us understand the weak dependence of the low energy excitations on the con- 
centration X. The large prefactors for J4 and Jg indicate that modest changes in J4 and Jg 
are sufficient to offset the large changes in Ji found within our model. At the X point, Eq. |^ 
reduces to 

huj{q^ = ^) - + 8-^2 + 2Jg. (5) 

As the dependence of J2 on x is weak, the changes at the zone boundary are driven by Ji 
and Jg. Unlike in the low q regime, the prefactors of Ji and Jg are equal in this case. Since 
the decrease in Ji is much larger than the increase in Jg, the energy at the X point decreases 
as the hole concentration is increased. Another useful information about the flattening of 
the dispersion beyond half way to the zone boundary is given by comparing Eq. 5 with 

= ^) ~ J1+4J2 + 2J4 + Jg. (6) 

The energies at q^ = ^ and qx = are comparable when Ji ~ 2J4. 

We constructed simpler models to make quantitative estimates for the contribution from 
the t2g electrons and that from the Cg electrons. Model B(C) includes eg {t2g) orbitals on 
the Mn atoms and all p orbitals on the oxygens. As the density of states (DOS) obtained 
within model A suggests partial occupancy of the minority spin t2g bands with ~ 0.175 
electrons even for the undoped case, the Fermi energy {E-p) of model C in the undoped case 
was adjusted so that the minority spin t2g bands had 0.175 electrons. As a consequence, 
the number of holes in the majority spin tg states of model B in the undoped case should 
be larger by 0.175 than the case when the minority spin t2g states are not occupied. The 
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d partial DOS for models B and C along with the result for model A which considers all 
d orbitals on the Mn atom are shown in Fig. 3. The reduced models (B and C) are found 
to give a good description of the respective d partial DOS of Cg and t2g symmetry within 
model A. Further justification for the treatment of the contributions from Cg and t2g states 
separately is given by the fact that the dominant contributions, (Jj, i=l,4,8), are all for 
the pairs along < 100 > for which there is no mixing of the two states in the exchange 
coupling. The spin wave dispersion was calculated in the reduced models separately. The 
dispersion along the FX direction for model B is shown in the panel b (inset of Fig. 3a). 
y denotes the number of doped eg holes with reference to the half filled majority spin eg 
band. (Equivalently, 1 — y is the number of electrons in the eg band.) x in the parentheses 
indicates the hole concentration in model A being equivalent to doping of divalent atoms. By 
considering the above situation, 7/=0.175 corresponds to the case of undoped LaMnOs, i.e., 
X — 0. As the number of holes increases from y = 0.175, the spin- wave energy at the X point 
steadily decreases. In model C, on the other hand, the dominant contribution from the t2g 
states to the exchange coupling is antiferromagnetic superexchange. The negative spin-wave 
energy for all q (Fig. 3d) is consistent with this expectation. However, small occupation of 
the minority spin t2g states produces a ferromagnetic double exchange contribution, z in 
Fig. 3d denotes the number of electrons in the minority spin t2g states. Clearly, doping of 
divalent elements reduces z so that the double exchange contribution diminishes rapidly as 
is clearly seen in the z dependence of the spin-wave dispersion in Fig. 3d. Here again the 
corresponding value of x in model A are indicated in brackets. The contribution from the Cg 
states (Fig. 3b) and that from the t2g states (Fig. 3d) for the common x value are added and 
the resultant spin wave dispersion shown in Fig.3e agrees very well with the one in Fig. lb 
(model A). This analysis suggests that the main source of the zone-boundary softening of 
the spin wave dispersion by doping of divalent atoms for x < 0.3 is the reduction in the 
ferromagnetic double exchange of the t2g electrons. On the other hand, in the doping range 
of a; > 0.3, the t2g states may simply act as a source of antiferromagnetic superexchange and 
further softening and flattening of the spin wave dispersion comes from the eg states. 
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The doped holes within our model have considerable oxygen p character, and the earlier 
results suggested that the itinerant oxygen band could modify the various exchange 
interaction strengths. It was pointed out that the role of oxygen p states in the superexchange 
interaction is not only to mediate the d — d transfer but also to make a direct additional 



contribution [0,^. However in this treatment the banding effect of oxygen p states was 
neglected. This assumption will not be justified for quantitative arguments if the p band 
width is comparable to the p—d energy separation, which is the case in our systems. In order 
to obtain information about the role of the oxygen p band, we made further simplification in 
the model B that the hopping between oxygen atoms was neglected, i.e., pp<7=pp7r=0 (model 
D). As the neglect of the hopping between the oxygen atoms could reduce the Cg bandwidth, 
the spin wave dispersions were calculated for two values of pda - the value (—2.02 eV) 
estimated already by the fitting (results shown in Fig. 4a) and an increased value of —2.25 eV 
(results shown in Fig. 4b). The dispersions shown in Figs. 4a and 4b are qualitatively similar 
to each other. In both cases, the spin wave energy at the X point increases with doping 
being in contradiction to the behavior observed in Fig. 1 and Fig. 3b. Since this model does 
not take into account the antiferromagnetic superexchange contributions coming from the 

degrees of freedom and affecting primarily the nearest-neighbor magnetic interactions, 
the ferromagnetic coupling Ji remains to be the strongest interaction in the system. In such 
a situation, the form of the spin-wave dispersion is close to the cosine-like. 

To analyze further the role played by the inter-oxygen hopping, the exchange interactions 
{Ji} were obtained for the cases corresponding to Fig. 3b and Fig. 4. The results are shown 
in Fig. 5. As is expected, the behavior of Ji is not affected so much by the p — p hopping. 
By comparing the results of model A in Fig.2 with those of model B in Fig. 5, we see that 
the doping dependence of J4 and Jg comes primarily from the Cg electrons. Neglect of the 
p — p hopping (model D) strongly suppresses J4 and Jg for x > 0.2, and zone boundary 
softening of the spin-wave dispersion become less pronounced. On the other hand, increase 
of J2 with hole doping is enhanced by neglecting the p — p hopping. The energy at X point 
increases in Fig. 4 (model D) with hole doping, because the variation in Ji in this case is 
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not enough to offset tfie sfiarp increase in J2. 

In summary, we fiave analyzed tfie low-temperature spin dynamics of the doped mangan- 
ites with tight-binding models. Our results provide some insight into the roles of the t2g and 
the Cg states and also suggest that the channel of hopping between the oxygen atoms strongly 
modifies the exchange interactions. Thus for the correct quantitative and sometimes even 
qualitative description, the simplifications made by model Hamiltonian approaches which 
consider only the Cg orbitals are questionable. 

We thank Prof. D.D. Sarma for useful discussions. Part of the programmes used here 
were developed in Prof. Sarma's group. The present work is partly supported by NEDO. 
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I. FIGURE CAPTIONS 



Fig.l The spin wave dispersions obtained by (a) LMTO calculations and (b) the tight 
binding approach (model A) along the symmetry directions FX, XM and MR shown as a 
function of doping. 

Fig. 2 The doping dependence of the exchange couplings Ji, J2, J4 and Jg between atoms 
at (o 0), (a a 0), (2o 0) and (3a 0), where a is the lattice parameter. 

Fig. 3 The (a) minority-spin and (c) majority-spin d partial density of states within 
models A, B and C. The spin wave dispersions along FX as a function of doping within 
models (b) B and (d) C are shown along with (e) the combined contributions of models B 
and C. y refers to the hole concentration in the majority-spin Cg band with reference to its 
half-filled case, z is the electron concentration in the minority-spin t2g band, x is the net 
concentration of the doped holes and is given hy x = y — z. 

Fig. 4 The dependence of the spin wave energies on the Cg hole doping y along FX within 
model D. The hopping between oxygen atoms and the t2g orbitals on the Mn atom have 
been left out of the model, (a) pda=-2.02 eV and (b) pda=-2.25 eV. 

Fig. 5 The variation of the exchange couplings Ji, J2, J4 and Jg with the Cg hole doping 
y. Open circles are for the case (model B) including the hopping between oxygen atoms 
and pda=-2.02 eV. Open and filled squares are for the cases (model D) without the hopping 
between oxygen atoms: pda=-2.02 eV (open squares) and pda=-2.25 eV (filled squares). 
The t2g orbitals on the Mn atom have been left out of the basis set. 
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